Rotation of a Bose-Einstein Condensate held under a toroidal trap 
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The aim of this paper is to perform a numerical and analytical study of a rotatin g B ose Einstein 
condensate placed in a harmonic plus Gaussian trap, following the experiments of [25|]. The rota- 
tional frequency f2 has to stay below the trapping frequency of the harmonic potential and we find 
that the condensate has an annular shape containing a triangular vortex lattice. As Q approaches 
u>, the width of the condensate and the circulation inside the central hole get large. We are able 
to provide analytical estimates of the size of the condensate and the circulation both in the low- 
est Landau level limit and the Thomas- Fermi limit, providing an analysis that is consistent with 
experiment. 
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I. INTRODUCTION 

The investigation of rotating gases or liquids is a cen- 
tral issue in the theory of superfluidity since they give rise 
to quantized vortices [H, E|- During recent years, several 
experiments using rotating atomic Bose Einstein con- 
densates have led to the observation of vortices. These 
condensates are usually confined in a harmonic poten- 
tial with cylindrical symmetry around the rotation axis 
z. Two limiting regimes occur depending on the ratio 
of the rotation frequency fl and the trap frequency to in 
the xy plane. When f2 is notably smaller than cj, only 
one or a few vortices are present at equilibrium 0,0]. A 
Thomas Fermi analysis can be performed to analyze this 
regime because the coupling constant describing the in- 
teractions is often large in the experiments [E When 
n approaches ui, since the centrifugal force nearly bal- 
ances the trapping force, the radius of the rotating gas 
increases, and the vortices arrange themselves on a tri- 
angular lattice @, IE IE E3] ■ A new class of phenomena in 
this regime of fast rotation is predicted in relation with 
Quantum Hall physics [E EI EE EE Ei EE EE E3 ■ In- 
deed the one-body Hamiltonian written in the rotating 
frame is similar to that of a charged particle in a uniform 
magnetic field and one can use the Landau levels struc- 
ture to analyze the ground state of the condensate and 
describe the properties of the lattice. 

In order to analyze the regime of fast rotation, one 
approach consists in adding a quartic potential to the 
harmonic potential. For this type of potential, the trap- 
ping force is always greater than the centrifugal force so 
that the regime > u> can be explored. The condensate 
is then seen to exhibit a more complex structure with 
regards to its density distribution and the arrangement 
of vortices [IE El, El EE EH, HE EE 13 ■ In particular, 
a multiply quantized vortex, or giant vortex, appears for 
large values of the rotational frequency Q and the con- 
densate is located within a thin annulus 2l| . When O is 
decreased from this situation, a circle of vortices exists 



inside the condensate [22 1. 

A number of experiments have been performed in 
which a laser beam is shone into an otherwise harmoni- 
cally trapped condensate 25, 5E 21, 3> thus creating a 
trapping potential of the form harmonic plus Gaussian. 
Often in experiments, the laser beam is weak, hence the 
Gaussian term is small and for the purpose of analysis 
can be expanded so that the resulting potential can be 
approximated by a harmonic plus quartic potential. A 
different approach to analyze these experiments is to con- 
sider the full harmonic plus Gaussian trapping potential. 

The aim of this paper is to perform a numerical and 
analytical study of a rotating condensate placed in a har- 
monic plus Gaussian trap as in the experiment [2HI . [3 • 
The specific feature of the Gaussian potential with re- 
spect to the quartic one is that the rotation frequency f2 
cannot get arbitrarily large but stays below u, the trap- 
ping frequency of the harmonic potential. We will show 
that according to the parameters of the system, the con- 
densate can either be a disk or an annulus. Furthermore 
we will show that as approaches u, the condensate 
always expands to become a large annulus with a vor- 
tex lattice inside the condensate and a large circulation 
within the central hole. This is in contrast to the har- 
monic plus quartic trap, which develops a giant vortex 
and a thin annulus. Using the Lowest Landau level (LLL) 
states, we will give an analytical description of the phe- 
nomena seen numerically. We estimate the radii R±, R2 
of the condensate and the circulation around the inside 
hole, of order R1R2, thus much bigger than that given 
by a uniform lattice (which would be R\). 

This paper is organised as follows. Section II contains a 
brief formalisation of the problem, introducing the energy 
functional followed by various numerical observations in 
Sect. III. The lowest Landau level analysis for the regime 
n close to u> is presented in section IV which provides the 
main analytical results of the paper. Finally, section V 
is devoted to extra computations in the Thomas-Fermi 
regime. 
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II. FORMULATION 

A two-dimensional Bosc-Einstcin condensate trapped 
at absolute zero temperature can be described by a 
macroscopic condensate wave function (order parameter) 
The ground state of the rotating system is determined 
by minimizing the energy functional E' = E — £IL Z where 
L z =$*[z'(rx p)]^ is the z component of angular mo- 
mentum along the rotation axis (for linear momentum 
p). The energy functional, in the frame rotating with 
angular velocity O is then 



2m 1 



+ ^.\^-nL z dv, 
(i) 

with r 2 = x 2 + y 2 and where the integral is carried out 
over the spatial domain V. The trapping potential is 
composed of a harmonic plus Gaussian term 

V tr = Vq exp (-2r 2 /w 2 ) + ^mu 2 r 2 . (2) 

When the atoms are assumed to occupy the ground state 
of the harmonic oscillator in the z direction, with en- 
ergy huj z /2 and extension a z = y/ h/muj z , suppression 
of the condensate in the z direction is allowed provided 
the characteristic energy Huj z is very large in comparison 
with the other energy scales. Here uj z is the frequency of 
the confinement in the z direction. The two-dimensional 
coupling parameter is then Uq = y/8nh 2 a s N/ma z for N 
identical atoms with s-wave scatting length a s [13] • 

The system can be nondimensionalised by choosing 
u>, huj and ^Jh/{mui) as units of frequency, energy and 
length respectively. Thus, on defining a non-dimensional 
coupling parameter g = mUo/h 2 , the energy functional 
takes the non-dimcnsionaliscd form 



E' 



i|V?A| 2 + V(r)H 2 + |H 4 -^ 2 



for external toroidal potential trap 



V(r) = Ae 



-l 2 r 2 1 2 

+ 2 r ' 



dV, (3) 



(4) 



with A = Vo/fkv and inverse waist I — (2^/mwwg) 1 ' 2 . 
The energy functional (J3j) is subject to the normalization 



\ip\ 2 rdV = 1. 



(5) 



In this scaling, large rotation implies that fl gets close 
to 1. Note that in experiments I is often small so that 
the potential V(r) in Eq. §Z§ can be expanded to give 



V(r)~ ^(l-2Al 2 )r 2 + ^Al 4 r 4 , 



(6) 



served |21[ . However in this paper we retain the toroidal 
potential given by Eq. ((4]) for the numerical and analyt- 
ical analysis. 

We will perform a full numerical analysis of the exper- 
imental case of [25|] , which will lead us to a numerical 
and analytical description of several model cases which 
prove to be different from the harmonic plus quartic trap 
considered in 21, [^j. In particular, as gets close to 



1, the condensate has an annular shape, its width always 
becomes large and a vortex lattice is present with a circu- 
lation inside the annulus. We are able to estimate these 
various quantities. 



III. NUMERICAL OBSERVATIONS 
A. The Effective Potential 

When the condensate is put into rotation, the effective 
trapping potential to be considered is not given by Eq. 
but is instead given by 



V, 



eff 



V(r) 



(I) 



Therefore according to the values of A, I and fl, this 
effective potential can produce either a disk condensate 
or an annular condensate. To see this, notice first that 
the effective potential has a minimum that occurs for 
r = 7'o > given by 



provided 



^ log v 



2A1 2 
I -CI 2 



2Af 



ft 2 



> 1. 



(8) 



(9) 



from which a critical frequency around 1 — 2A1 2 is ob- 



If q > 1 , the effective potential has a local minimum and 
it can lead to two different situations, either the conden- 
sate is a disk or an annulus. For existence of an inner 
boundary, we must have q > 1 + S for some positive (not 
necessarily small) 8. As ft — > 1, q > 1 + 8 is always 
satisfied and so that an inner boundary is created. The 
determination of the value of 8 is not readily obtained 
as it depends on the normalisation condition ([5]). Con- 
versely if q < 1, then the condensate is always a disk. 

Figure [1] shows three examples of the effective poten- 
tial (O plotted against radial distance from the cen- 
tre of the condensate along constant 9 for the param- 
eters {g,A,l} = {100,25,0.03}, {1000,10,0.75} and 
{500,100,0.9} with = 0. In the first parameter set 
q < 1, the condensate is a disk and the density maximum 
is at the centre. In the second parameter set, there is a 
local density minimum at the centre of the condensate, 
but the condensate is still a disk. The third parameter 
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set displays an inner boundary and the condensate is thus 
annular. 
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FIG. 1: The effective potential given by Eq. J7} as a function of 
radial position is shown in the left column. The associated density 
(p = |V>| 2 ) plots taken along y = are shown in the right column. 
Three parameter sets considered are (a) {g,A,l} = {100,25,0.03}, 
(b) {1000,10,0.75} and (c) {500,100,0.9} all with Q. = 0. Dis- 
tances are measured in units of \JKj (raw), density in units of muj/h 
and potentials in units of fLo. 



The effective potential plotted in Fig. [T] only consid- 
ers a non-rotating condensate, SI — 0. When the con- 
densate is placed under rotation, vortices form and the 
shape and size of the condensate alter. Numerical simu- 
lations on the Gross-Pitaevskii equation are carried out 
to explore the effect of fi on a range of parameter sets 
for {<7, A, I}. The two-dimensional, dimcnsionless, Gross- 
Pitaevskii equation comes directly from the energy func- 
tional ((3]) and is 



dip 



1 



dip 



dx 



for V(r) given by (g]). 



dip'' 

dy, 
(10) 

The Gross-Pitaevskii equation is 

I) by 



solved numerically in imaginary time (see [2 ll 
evolving an initial wavefunction for a range of values of 
SI < 1 in order to find the ground state. Three cases 
of interest, which summarise the numerical results well, 
are presented below. The three reported parameter sets 
are; {g,A,l} = {955.95,24.83,0.07}, {14,1000,5} and 
{500,60,0.1}. 



B. The Experimental Case of Bretin et al. [2 

A natural case to numerically simulate is that consid- 
ered experimentally by Bretin et al. pjj where a har- 
monically trapped condensate is put in rotation with a 
weak laser beam shone at the origin, modeled by a Gaus- 
sian term. This experimental case can be described by 
a two dimensional system as explained in the introduc- 
tion, using the dimensional reduction which leads to the 
definition of Uq. The experimental values of [25| cor- 
respond to {g, A,l} = {955.95,24.83,0.07}. A series of 
contour plots of the density are shown in Fig. [2] (see Fig. 
1 of [25| , with the appropriate rescaling of rotational ve- 
locities, so that the fl^v = 60 of [25| corresponds to 

SI = 0.795 in this paper and Q st { r = 69 corresponds to 
SI = 0.914. The rotational velocity is calculated from 
[25| using the value of the frequency in the x and y direc- 
tion u>± /2tt = 75.5-ffz and not the second stirring phase 
frequency ui±/2ir = 64.8iJz). 

For the slow rotational velocities (0 < 0.5) of Fig. [21 
the condensate is a disk with a small number of vortices 
(12 vortices are present when SI = 0.5; see Fig. &b)) 
with the vortices forming a triangular lattice. As the 
rotational velocity is increased, the radius of the outer 
boundary increases while more vortices are accommo- 
dated into the condensate. The dynamics here mimic 
those observed in harmonic traps. 

Above some angular velocity the density at the cen- 
tre of the condensate begins to deplete and eventually an 
inner boundary, hence an annulus, is created. In the ex- 
periments of [25j, a density minimum at the centre first 



occurred for SI ~ 0.874 (corresponding to SI 



(2) 



66; 



see plate (c) of Fig. 1 in (25[). The numerical simula- 
tions here suggest the onset of the density minimum to 
be SI ~ 0.887; see plate (d) of Fig. [2] where the density 
minimum first appears. Clearer pictures of the develop- 
ment of the depletion of density at the centre can be seen 
in plates (e) and (f). It distorts the vortex lattice in much 
the same way that the outer boundary does. In Sect IV 
we note that, under the LLL approximation, the vortex 
lattice inside the central hole is distorted from a regular 
vortex lattice such that the number of zeros in the hole 
is given by R1R2, to get a number of order R\, where 
i?i,i?2 are the inner, outer radii. Increasing the angu- 
lar velocity still further, thus exploring the fast rotation 
regime f2 — ► 1, details how the density at the centre of 
the condensate continues to diminish until for SI ~ 0.92 
a central hole develops (see plate (g)) and the conden- 
sate becomes annular. The central hole grows rapidly; 
for SI = 0.93, the central hole is large and there is a cir- 
culation equivalent to 11 vortices (see plate (h)). Our 
simulations have been carried out up to ^ 0.95. These 
higher angular velocity simulations suggest that both the 
outer and inner radii and also the width of the condensate 
increase in size as O — > 1. 
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in the vortex lattice and requiring a longer time of con- 
vergence for the numerical simulations, must be one of 
the reasons that explain the experimental difficulties in 
observing the condensate at these rotation frequencies. 



C. The Annulus 

Manipulating the values of the parameters {g, A, 1} can 
have the effect of altering the shape of the condensate. 
Here we will consider a parameter set that creates an 
annular condensate that is present for all CI < 1. The pa- 
rameter set is thus chosen to be {g, A, 1} — {14, 1000, 5}. 
A selection of contour plots for various angular velocities 
are given in Fig. [3l The choice of this parameter set is 
to illustrate the behaviour as f2 gets close to 1. 
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FIG. 2: (Color online) Density profiles of a rotating condensate 
with {g,A, 1} = {955.95,24.83,0.07} corresponding to the (non- 
dimensionalised) experimental values of [25ll for (a) Q = 0.25, (b) 
fi = 0.5, (c) n = 0.874, (d) Q = 0.887, (e) Q = 0.901, (f) Q = 0.914, 
(g) f2 = 0.92 and (h) Q = 0.93. Distances are measured in units of 
\Jh/ (mu) . 

The experiments of Bretin et al. [lH provide an ex- 
ample of the transition from a disk condensate with the 
density maximum at the center to a disk condensate with 
a local minimum at the center. It is reasonably safe to 
assume that if the angular velocity could be further in- 
creased in the experiments, the condensate would become 
annular, with a large persistent current. The depletion of 
density, which occurs for Q, > 0.887, creating a distortion 



FIG. 3: (Color online) Density profiles of a rotating condensate 
with {g,A,l} = {14,1000,5} for (a) U = 0.25, (b) U = 0.9, (c) 
f2 = 0.99 and (d) Q = 0.994. Distances are measured in units of 
yh/ (maj). 

For low rotational velocities, see Fig. [3^a), the conden- 
sate does not contain vortices in the annulus. However 
the closer Q approaches unity, both the radius of the in- 
ner and outer boundaries increase, but so too does the 
width of the condensate. When Q = 0.9 (Fig. G^b)), the 
condensate still does not contain any vortices, but for 
il = 0.99 (Fig. [UJc)), the condensate contains two com- 
plete rings of vortices. For fl = 0.99, there is a multiply 
quantised vortex at the centre of the condensate provid- 
ing a persistent flow with a quantum of circulation v = 3. 
The phase profiles also show that there are further singu- 
larities of phase ('invisible' vortices) in the outer regions 
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of the condensate. 

Note that, for all O, the condensate is always an annu- 
lus with the inner and outer radii both increasing as f2 
increases. Furthermore the width of the condensate also 
increases so that a thin annulus is never created. The 
increase in size of the condensate for D, > 0.9 is marked. 
This parameter set explicitly shows the presence, at large 
il, of a central hole containing circulation together with 
a vortex lattice in the bulk of the condensate. 

The choice of this parameter set, especially the value 
of g, is specifically chosen with reference to the lowest 
Landau level (LLL) analysis in Sect. IV. As will noted in 
Sect. IV, to use the LLL approximation requires g(l — 
il 2 ) to be small. As a consequence, with an eye on the 
capability of the numerical simulations to resolve at f2 
close to 1, it necessarily forces g to be not too large, 
though the main features are preserved while increasing 
g. We note here that an annular condensate, existing at 
all fi, can be created for a wide range of values of g. 
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D. Density Dip at the Center 

As a final numerical example, one can consider the 
parameter set {g,A, 1} = {500,60,0.1}. The choice of 
these parameters actually forces, for = 0, the ground 
state to have a local non-zero minimum of density at the 
centre. For these parameters the density maximum is 
located at r = 4.2 when = and the condensate is a 
disk; a contour plot of the condensate at fl = is given 
in Fig. [4] along with a selection of other plots for various 
angular velocities. 

Vortices first appear close to the maximum of density 
instead of close to the center of the condensate. As il is 
increased, the vortex lattice develops close to this initial 
circle of vortices. The increase in size of the condensate 
is visible, as is the depiction of density at the centre of 
the condensate, which turns into a hole: an annulus is 
formed. 



IV. Lowest Landau Level Analysis 

We turn to the analysis of the ground state of the en- 
ergy ([3]), where V(r) is given by Eq. (j4j and fl tends to 1. 



We will use a Lowest Landau level analysis 
Recall that the spectrum of the Hamiltonian 

H Q = -iv 2 + ir 2 -r>L 2 , 



11, 17, 3 



(11) 



has a Landau level structure. The lowest Landau level is 
defined as 

f{x + iy)e-%( x2+y2 ), / analytic. (12) 



FIG. 4: (Color online) Density profiles of a slow rotating conden- 
sate with {g,A,l} = {500,60,0.1} and (a) Q = 0, (b) O = 0.225, 
(c) f2 = 0.35 and (d) Q = 0.5. Distances are measured in units of 
y/h/ (mid). 



For such functions, < Hqi/}, ip > can be simplified (see 
[HI) so that E' = Q + E LLL (ip) where 



E 



LLL 



WO 



1 



V{r) - £fiV I |V^| 2 + J WdV. (13) 



The minimization of Eq. (|13| without the analytic con- 
straint provides the Thomas-Fermi profile for the coarse 
grain density: 



h/f^PTF := A - 1 (V(r) - l^V 



(14) 



Such a function can be recovered in the LLL by the pres- 
ence of the vortex lattice and this only changes the co- 
efficient g into bg where b is the Abrikosov parameter as 
we explain further below. 

We recall that the orthogonal projection for a general 
function -0 onto the LLL is explicit [2^, [3^1 and given by: 



n 



LLL 



e -H\ z \ 2 - 2z z' + \z'\ 2 )^( X ',yl)dx'dy', 



(15) 



where z = x + iy and z' = x' + iy' . If an LLL function ip 
(i.e ip satisfies Eq. (fT2)l ) minimizes the energy Eq. (| 13[) . 
it is a solution of the projected Gross-Pitaevskii equation: 



n 



LLL 



[(V(r) - ftV/2 + 5 |Vf -n)i/>]= 0, (16) 
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where p is the chemical potential. 



A. Disc Condensate 



When is close to 1, and V(r) is the harmonic 
plus Gaussian potential, Eq. (|16|) can be approximated 
by TIlll [<?|V'| 2 V , ] = LLil>, which is the equation of the 
Abrikosov problem (see [11, 31, 32|). A solution can then 
be constructed using the Theta function (see [32| for the 
details): 



{x,y;r 



) = e M* 2 -N 2 )e 



(17) 



where r = tr+iti is the lattice parameter. The zeroes of 
the function <fi lie on the lattice . — (Z © It) and \<$>\ is 



periodic. The optimal lattice, that is the one minimizing 
p{r) = J \4>\ 4 /(J \4>\ 2 ) 2 , is hexagonal, which corresponds 
to r = e 2l7r / 3 (the integrals are taken on one period). For 
r = e 2 "^ 3 , (j,(t) = b~ 1.16. 



As in [32], we can construct an approximate ground 
state by multiplying the solution, Eq. (fT7|) . of the 
Abrikosov problem by a profile p varying on the same 
scale as /?tf, which is large. Since this product is not 
in the LLL, we project it onto the LLL and define 

v = Rlll 



ergy of v yields 



V(r) 



Estimating the en- 



gb 



OV) p(x,y) + ^-p(x,yf 



This computation assumes that <f> and p do not vary on 
the same scale, hence the integrals can be decoupled. 
Then, minimizing with respect to p implies that p must 
be a Thomas- Fermi profile (Eq. (fT4]0 with g changed 
into bg: 



p{ x , y ) = \-L(v{r)- l -tfr* 



(18) 



This approximation is valid provided the energy ob- 
tained, £'lll(w), is much smaller than the gap between 
two Landau levels which is of the order unity. 



In the case of the harmonic potential V(r) = r 2 /2, 
the Thomas Fermi profile provides a disk condensate of 
radius R = [Agb/{ir{\ — fi 2 ))] 1 / 4 which is large when Q 
gets close to 1. Moreover, -Elll(v) is of order sj g(l — fl 2 ) 
which is indeed small when O is close to 1 and g is not 
too large, so that the LLL approximation is satisfied. In 
the case of the toroidal potential ([4]) , we have to compute 
the Thomas- Fermi profile from Eq. (fT8|) and discriminate 
whether it is a disk or an annulus. Then we have to 
check whether the LLL approximation is justified, that 
is whether -Elll(v) is small. 



Suppose that the condensate is a disk (recall that this 
requires that q is not too large), so that only an outer 
boundary exists and write 



gbM 2 = p + -(n 2 -l)r 2 



Ae~ 



(19) 



To find approximations to the radius of the outer bound- 
ary one can proceed by taking Eq. (|19|) and the normal- 
ization condition (J5J) which provide the necessary starting 
equations in order to compute p. A check on the validity 
of the LLL approximation is to verify, in the limit £1 — > 1, 
that the chemical potential p is small. 

To begin, substitute the density, \tp\ 2 , from Eq. (fl9|) 
into the normalization condition ([5|) and integrate over 
the domain between and i?2, where i?2 is defined as the 
radius of the outer boundary. It follows that 



gb 



pR\ 



1 



A 



(n 2 -l)Rt + -(e 



1), (20) 



which explicitly contains the chemical potential p. To 
remove p from the calculations, one can note from Eq. 
(HE) that 



p-V eff (r) 



from which 



gb 



p = Ae 



-(tt 2 



Using Eq. (J22J) in Eq. J2QJ) gives 



^ = hl-n 2 )Rl + A 

7T 4 



Rle- lR '- 



l)R 2 . 



(21) 



(22) 



l Rn 



1) 



(23) 



At this stage we introduce the parameter p, defined as 

gi* 



p 



ir(l-n 2 



(24) 



We are interested in the limit fl close to 1, i.e. p 3> 1, 
which corresponds to l 2 R 2 ^> 1. Then the exponential 
terms in Eq. (|23[) immediately disappear and it follows 
that 



i?2 



(1-tt 2 ) [n 



gb 



1/4 



and 



(l-Q 2 ) 



gb 



A 
I 2 



1/2 



(25) 



(26) 



This implies that I R% is large and p is small since fi is 
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close to 1. Recall that the assumption of a disk conden- 
sate requires q to be bounded, q < 1 + 5, and thus (as i5 
is not large) 



Al 2 < 



(l + <5)(l-ft 2 ) 



(27) 



so that the ratio Al 2 has to be small in this regime. 



B. Annular Condensate 



Expressions (|25|) - ([26|) arc valid strictly when there is no 
inner boundary which occurs when q < 1 + 5. 1£ q > 1 + 5 
then an inner boundary, at r — Ri, develops. 

One again starts from Eq. (TlT)|) , but this time the in- 
tegration is taken over the domain between R\ and R 2 . 
It then follows that 



% = ^Rl-R\)+\itf-l){Rt-R\)+^ (e-^ 



(28) 

which again explicitly contains the chemical potential p. 
To remove /x from the calculations, we note from Eq. (|19p 
that 



V-Veff{r) 



fjb 



\R1.R2 



0, 



from which we recover Eq. 



and 



-l Rn 



D -l 2 Rj _ \ ic 



= —{W-l){Ri-R{). 



(29) 



(30) 



and Eq. (JSHJ) into Eq. 
and writing the 'area' of the annulus as X = R 2 — R\, 



Upon substitution of Eq 
and w: 
we get 



gb = ^X(ft 2 - 1) 



Xe 



-i A x 



-PX 



1 



1 



X 

T 



(31) 



Since p is large (because ft is close to 1), we can assume 
that l 2 X is large so that the first term in the squared 
brackets of Eq. (j3Tj) is negligible, which leaves 



Thus one obtains an expression for X, 



X 



Agb 



7r(l - ft 2 )' 



(34) 



so that l 2 X is indeed large when ft gets close to 1. Fur- 
thermore, from ([3T)]) 



-l 2 R-l _ 



(ft 2 - 1)X 



2A(e 



1)' 



(35) 



which, on using the derived approximation for X, Eq. 
([34)l . implies that 



Ri 



In A 



gb(l - ft 2 ; 



-,1/2 



(36) 



giving an expression for the radius of the inner bound- 
ary. The value of the chemical potential can then easily 
be found by taking Eq. (|2"9"|) evaluated at with the 
expression for Ri following directly from Eq. (|36l) ; 



1 



(l-ft 2 )i? 2 + Ae 



-VRi 



gb(l - ft 2 ) 



provided 



In A 



gb(i-n 2 ) 



Agbl 4 



?r(l "ft 2 ) 



1/4 



(37) 



(38) 



In the limit ft — > 1, we see indeed that n — > 0, thus 
justifying the LLL approximation provided g is not too 
large. The radius of the outer boundary can be found by 
taking Eq. (|22J) with fj, given by ([37]). Then 



R 2 



4gb 



tt(1 - ft 2 



1/4 



(39) 



which is just y/X, implicitly implying through ([38| that 
R 2 3> R 2 - 

To find the width of the condensate, d = R 2 — R\, it is 
necessary that i?i be neglected in front of R2, i.e. that 
R1/R2 *C 1 which again is equivalent to Eq. ([38]) being 
satisfied. Thus, the width of the condensate d is found 
from 



gb= ^7rX 2 (ft 2 - 1) 



1 

PX 



(32) 



Now, as l 2 X is taken to be large, one can neglect l/l 2 X 
in front of the factor of a half in the square bracket. 
Formally, this corresponds to the following being satisfied 



1 7T 

^ft 2 > ~gP' 



(33) 



X = 
^d = 



{R 2 -R 1 ){R 1 +R 2 ) 

X X _( 2gb 

R1+R2 ~R~2~ \k{1-S1 2 ) 



1/4 



(40) 



Notice that both the width of the condensate d and the 
inner R\ and outer boundaries R 2 get larger as ft in- 
creases and eventually tend to infinity as ft — ► 1. How- 
ever, while the width of the condensate and the outer 
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TABLE I: A comparison, for parameters {g, A, 1} — 
{14, 1000, 5}, of the values of the inner radius Ri, outer radius 
i?2 and quantum of circulation v for fast rotation calculated in 
the LLL and given by Eq.'s (J36J) , ^39j and (|42j) respectively. 
Numerical values are provided (subscripts n) as a compari- 
son. The chemical potential fi is also given, calculated from 
Eq. (£7J. 



n 


Ri 


Rln 


i?2 


Rln 


V 






0.9 


0.53 


0.66 


3.23 


3.25 


1.70 


2 


0.99 


0.99 


0.57 


0.69 


5.68 


5.71 


3.22 


3 


0.32 


0.994 


0.58 


0.70 


6.45 


6.38 


3.72 


3 


0.25 



boundary grow like (i - n 2 )- 1 / 4 , the inner boundary 
grows like ln[(l - ft 2 )" 1 / 2 ] 1 / 2 . Thus, in the limit fi -> 1, 
the condensate forms an infinitely thick annulus with 
both the inner and outer boundary tending to infinity 
and with an infinitely large central hole. 

The Cauchy formula allows us to compute the number 
of zeroes v{Ri) inside the disk of radius R\: 



V. Thomas-Fermi approximation 

In the case of the experiments of Bretin et al. [25|, 
is not very close to 1 and g is large, so the LLL analysis 
of Sect IV docs not adequately describe this experiment. 
In order to describe it better, and since g is large, we can 
use the Thomas-Fermi (TF) approximation. In the TF 
approximation, the vortex cores are a small perturbation 
with respect to the density profile and yield similar equa- 
tions to those previously derived in Sect IV, except that 
now there is no factor b entering the TF density profile 
of Eq. ([TO]) since the vortex cores are small. These com- 



putations are similar to those in 2l| , however due to the 
Gaussian trapping potential used here, the condensate 
becomes a thick annulus, with many vortices, instead of 
a thin annulus with no vortex as in 2l|. This changes 



significantly the approximations and requires the analy- 
sis of several cases according to the magnitude of l 2 R 2 , 
R being the radius of the condensate. 



A. Disk condensate 



2tt 



2tt 

d0 



gi — 

J e R ? z e~ 



(41) 

where (j> comes from (fl~7"l). Using a Laplace method to 
evaluate the integrals, we see that for large Ri, since 
R1/R2 tends to 0, 



v{R 1 )~R 1 R 2 . 



(42) 



Note that a regular lattice in a disk of radius R\ would 
give v(Ri) ~ R 2 , which is much smaller. 



If there is no inner boundary then the Thomas-Fermi 
approximation leads to 



IV 



(43) 



We substitute the density \ip\ 2 from Eq. (|43[) into the 
normalization condition ([5]) and integrate over the do- 
main between and R2 , where as before R2 is defined as 
the radius of the outer boundary. It follows that 



^=, l R 2 + 1 { n 2 -l)Rl + ^(e- 12 ^ 

7T 4 L z 



1), (44) 



which explicitly contains the chemical potential fi. To 
remove /i from the calculations, one can note from Eq. 
(1431) that 



C. Summary of the main results 



V>-V e ff(r) 



9 



M 2 L =0, 



(45) 



We have derived that as fl approaches 1, the conden- 
sate has an annular shape with a triangular vortex lattice 
inside. Eq.'s (36]), (f39j) , (37|) and fl42|) provide the inner, 
outer radii of the condensate, the chemical potential and 
the circulation in the inside hole. In the LLL approxima- 
tion, we must have y/ g(l — SI 2 ) <C 1. Therefore in the 
limit — >• 1, the parameter set {g, A,l} = {14, 1000,5}, 
described in Set. IIIC and with a series of contour plots 
(Fig. [3]) showing that the condensate is always annu- 
lar, satisfies this condition. Table [I] gives a comparison 
between the numerical and analytical values for this pa- 
rameter set and provide a good check on the estimates 
derived in Eq.'s (0SD, (07]), (|3"9 l and 



from which we get 

/.t = Ae7 

which is the same as 
yields 



5 ( n " 1)^2, 



(46) 

Using Eq. (gBJ) in Eq. flUD 



9_ 

IT 



1 



(1 -n 2 )R^ + A 



Rie 



(47) 

If I is taken to be small as in the experiments, we can 
assume that p is small and thus that l 2 R\ is small. Then 
the exponentials in l 2 R\ can be expanded and it readily 
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TABLE II: A comparison, for parameters {g,A,l} — 
{955.95, 24.83, 0.07}, for Q such that the condensate is a disk. 
The value of the radius of the condensate for different rota- 
tional velocities calculated asymptotically from Eq. ((48} un- 
der the assumption of the TF approximation and calculated 
numerically is shown. 



n 

o.i 

0.25 
0.5 
0.795 
0.821 



R2 (asymptotically) i? 2 (numerically) 



6.29 
6.40 
6.87 
8.97 
9.44 



6.18 
6.19 
6.70 
9.00 
9.49 



from the TF approximation of \tj)\ 2 given by ([43]) . but 
not exactly since the presence of the circle of vortices 
has an influence on its shape. A computation similar to 
that in [22| would be required to determine p in this set- 
ting. Once this is done, one should be able to use the 
results of (33|, chapter 3, which states that the radius 
where the circle of vortices appear is given by the loca- 
tion where the function ' p reaches its maximum, where 
£(r) = J sp(s) ds, R being the radius of the condensate. 
This follows from an expansion of the Gross Pitaevskii 
energy, where the leading order is given by the energy 
of vortices of order p(r) In £, where £ is the scattering 
length, minus the L z term, which can be estimated as 
—fl((r). Thus, the lowest f2 for nucleation of vortices is 
achieved when there is a radius where £/p is minimal. 



follows that 



jZ 4 fl|(l -q) + -ql 6 R 



2' 



(48) 



where we have used the expressions for q and p (Eq.'s ([9|) 
and (|24|) respectively) . Equation |48|) is a cubic equation 
in l 2 R 2 . When q < 1, the last term in Eq. (J48]) can be 
neglected and we get 



l^l-q) 7r((l - Q 2 ) - Al 2 ) ' 



(49) 



When q becomes larger than 1, then 1 — q < 0, but since 
we assume that we have a disk condensate, we recall that 
1 — q is small. In such situations, the last term in Eq. 
(|48|) cannot be neglected and is crucial to get the correct 
sign on the right hand side of |48|) . 

The chemical potential is derived from Eq. (|46|) 



A 1 



1 



1 



Rl(l-&)(l-q) + -Al' i Rt(l + q), 



(50) 

with the value of R2 obtained from Eq. (|48|) . This al- 
lows us to check that glV^O)! 2 ~ p — A is large, provided 
1 — SI 2 is not too small, g is large, I is small so that 
p is small. Hence, this justifies our use here of the TF 
approximation. This is in particular the case for the pa- 
rameters {g,A,l} = {955.95,24.83,0.07} corresponding 
to the experiments of [25j |. 

Table |TT] gives a comparison between the numerical val- 
ues and analytical values (calculated from Eq. (|48|) ) for 
the boundary of the condensate when the condensate is 
still a disk (f2 < 0.92), with the agreement found to be 
extremely good. 

A specific feature of our numerics in the case {g, A, 1} = 
{500, 60, 0.1} is to find that for certain values of the pa- 
rameters, since the maximum of the density is not at 
the origin, that vortices appear on a specific circle rather 
than at the origin (see Figure Sfb)). We call p(r) the 
average of \ip\ 2 on a circle of radius r. Then p is not far 



B. Annular condensate 

When the condensate is an annulus, then in the TF 
approximation, we need to take into account the quan- 
tum of circulation v in the inner hole of the condensate. 
We assume that ip = |-0|e Ily-S , where S is the phase. Note 
that this is not needed in the LLL approach, because the 
circulation is incorporated into the LLL wave function 
and we can compute it directly from (|41|) . 

The TF density expression (|43|) is adjusted to 



(n 2 - l)r 2 



Ae 



(51) 



where p, = p + Qv. The TF density (|5ip and the nor- 
malisation condition ([5]) provide the starting points from 
which approximations to the values of the radii of the 
condensate boundaries and thus the width of the con- 
densate can be found. It follows then that 



(Q 2 l)r 2 



-Ae 



(52) 



and eliminating p from Eq. (|52[) gives 



E>2 r?2 



(l_0 2 ) + ^ e -^(e- i2 



1) 



, (53) 



where X is the 'area' of the condensate and is again de- 
fined as X = R\ — R\. From integration of the normali 
sation condition j5]) between Ri and Ri it ensues that 



9_ 

IT 



1 



1 



■(1 - Q 2 )X(R( + R^ - -(1 - n^RfR 1 In 



R 2 2 

Ri 



Ae - l2R l(e- l2x -1) 
AXe- l2x e- l2R l 



1 

/2 



1 2 \n ' 2 



X 



Ri 



(54) 
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A further expression that connects the quantum of cir- 
culation and the inner and outer radii is readily obtained 
from the minimisation of the free energy per particle, 
F = E' — jLtiV, where N is the number of bosons in the 
condensate [2l|. The following integral identity results: 



g£l =2w - [g\^\ 2 ] dr. 



(55) 



Integration of Eq. (|5"B")) using Eq. (|43[) . together with the 
expressions for /2 and v above, results in 



TABLE III: A comparison, for parameters {g, A, 1} = 
{955.95, 24.83, 0.07}, for f2 such that the condensate is an an- 
nulus. The values of the inner radius Ri , outer radius R2 and 
quantum of circulation v are calculated by a TF analysis and 
are given by Eq.'s (|62p . (|59|l and (|61|l respectively. Numerical 
values are provided (subscripts n) as a comparison. 



n 


Ri 




i?2 


Rin 


V 


v n 


0.92 


1.00 


0.00 


12.94 


13.30 


1.30 





0.93 


3.72 


3.25 


13.52 


13.80 


1.37 


11 



gn 



(1 



tt 2 )(i?? 

He 



2 

A Rj -l 2 Rh -l 2 X 



R 2 ) 



X 

Ae~ l2R i(e- l2x 



Ae~ 



lnlS 



1 



2A 



Ri 

R 2 



x(i - n 2 



dr. 



Ri 



(56) 



In order to proceed it becomes important to estimate the 
last integral in Eq. ([56]) according to the size of l 2 r 2 . The 
following sections detail two possible limits. 

l 2 R\ small. If l 2 R\ is assumed to be small, then it 
follows that l 2 R\ is small as well. Thus expanding the 
exponential terms in Eq. (|53|) up to terms of order l^R% 
gives 



R 2 R 2 [{l-n 2 )(\-q)+Al i R 2 ] 



(57) 



where we have assumed that R\ ^ R\, which implies 
that we can write exp(— l 2 X) ~ exp(— l 2 R%). In the case 
of an annular condensate, we have that q > 1, so that 
the first term in Eq. (|57|) is negative. Therefore in order 
for the right hand side of Eq. (|57|) to be positive, it is 
required that 



/i?2> AP 



Q 2 



1 

1 - - 

q 



(58) 



Expansion of Eq. ([54")) to order Z 6 i?2 (again assuming 
that R 2 ^ Rf) results in an expression for R2, 



1 



(1 - q)l 4 R 



Ri 1 



ql 2 ( 



Rl 



(1 - 9) V 3 
2R 2 \n{R 2 ) 



R\ Yn{R 



l 4 Rt(l-q) ql 6 R 6 2 
4 + 6 : 



(59) 



where we assume that R\ » R\ \n(Rl/R 2 ) ~ R\ ln(i?2)- 
Notice that (|59|) is actually the same approximation for 
i?2 as calculated for the disk condensate in the TF ap- 
proximation (see Eq. (|48|1 ). Finally we can consider Eq. 
(|56[) . Expanding the exponential under the integral in 



powers of l 2 r 2 and simplifying (|56p gives 
|(l-^)i?l [ln(i?l)-2] 



2\n{R 2 )-l 2 R 2 (ln(i?2) + l) 



+ -l 4 Rt(HR 2 ) + - 



(60) 



ln(i?|) 



-(l-n 2 )R 2 2 (l-q) + 2A 



where we assume that ln(i?2/^i) ~ ^ 1- It fol- 

lows that 



2gVt 
7T \n{Rl 



[R 2 2 (l-n 2 )(l-q) + 4A\ 



(61) 



and thus the inner radius, R\, is computed from (|57|) to 
be 



^i~"5- [(l-« 2 )(l-9) + ^ 4 ^] 



-1/2 



(62) 



There are now three equations (|59|) . ([BTjl and ([62]) for R 2 , 
v and i?i that describe the annular condensate in the TF 
regime with l 2 r 2 small. These provide a comparison to 
the numerical values corresponding to the parameters of 
Bretin et al. [2j|. A comparison for the case of £1 = 0.92 
and SI = 0.93 between the numerics (c.f. Fig. [Hg,h)) and 
analytical estimates is given in Table Unl The TF analyt- 
ical expressions describing the radii of the inner and outer 
boundary are reasonnably good. However the quantum 
of circulation (Eq. (frJTj) ) is not providing a suitable esti- 
mate. Indeed v is still small and the inconsistency could 
be as a result of the difficulty in numerically discriminat- 
ing between vortices inside the condensate and inside the 
inner boundary. 

We note that the Thomas Fermi approximation is jus- 
tified because the maximum of (?|'0| 2 is larger than 1 (nu- 
merically around 4). Larger values of fl which are not 
reached by the experiments would be better described by 
a LLL regime. Nevertheless, in the TF approximation, 
we can still analyze an intermediate case. 

l 2 R\ small and l 2 R 2 large. If we take l 2 R 2 small and 
l 2 R\ large then this implies that l 2 X is also large. Taking 
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Eq.'s (|53|) . (|54|) and (|56|) as the starting point we note 
immediately that (|53[) simplifies to 



boundary 



i? 2 i? 2 (l-fj 2 ), 



(63) 



2tt ln(i?; 



(69) 



while Eq. 



simplifies to 



5 

7T 



1 



n 2 )RlRl\n(Rl) 



fi 2 )^ 2 . [i? 2 -2i? 2 ln(i? 2 )] -Ai? 2 1 



1 



again for R 2 given by Eq. (1671) . 

Such a case is not reached experimentally and is just 
before the LLL regime. One could perform similar com- 
putations assuming that l 2 R\ is also large. 



VII. CONCLUSION 



l 2 R\ 



-(1-fF 

4 V 



^2 - 7a> 



where we have assumed that 



R 2 > Ri, 



R\ » 2i? 2 ln(i? 2 ), 



the first condition following directly from the assumption 
that l 2 R\ is small and l 2 R\ is large. Furthermore Eq. 



(|56|) simplifies as 

qVL \l , 

-(1 

2 V 

A 



n 2 )R 2 



2AR{ 



21 



X 



A + 2Aln(lR 1 ) 



\n(R 2 ) - R 2 (l - O 2 ) 
A 



-i*m 



PR 2 ' 



Motivated by the experiments of [251 ]. we provide nu- 
(^nicrical and analytical computations which describe the 
properties of a condensate placed in a harmonic plus 
Gaussian trap. We have seen that, however close f2 
gets to the harmonic trapping frequency, the condensate 
becomes a large annulus containing a triangular vortex 
lattice, contrary to that seen for a condensate with a 
quadratic plus quartic term where the width of the con- 
densate decreases. We estimate the circulation in the 
central hole which is higher than that corresponding to 
a regular lattice in this region. Also in a Thomas Fermi 
approximation, when the rotational velocity f2 is not too 
f close to the harmonic trapping frequency, we estimate the 
radii of the condensate and the circulation in the inside 
hole, in a way which is consistent with numerics. 



\{l-n 2 )Rl [ln(i? 2 )-2] 
hl-n 2 )R 2 2 \n(Rl), 



(65) 



since 
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-e dr 



ln(LRi) 



— 1~ R 



2c 



2l 2 R 2 



(66) 



in the limits l 2 R\ small and l 2 R\ large. Thus we get an 
expression for the outer boundary, from Eq. (|64[) . as 



Ri 



(1 - O 2 ) \ir I 



A 



-l 1/4 



(67) 



and from Eq. (|65p we get an expression for the quantum 
of circulation 



2gVt 



ir(l-n 2 )R 2 ln(i? 2 )' 



(68) 



with i?2 found from Eq. (|67|) . Notice how the value of v 
gets large in the limit ft -> 1. Putting Eq.'s ([57]) and (|B"5)) 
into Eq. (|63[) we arrive at the expression for the inner 
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